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Abstract. Massive X-ray binaries are usually classified by the properties of the donor star in classical, supergiant and Be X- 
ray binaries, the main difference being the mass transfer mechanism between the two components. The massive X-ray binary 
4U 2206+54 does not fit in any of these groups, and deserves a detailed study to understand how the transfer of matter and the 
accretion on to the compact object take place. To this end we study an lUE spectrum of the donor and obtain a wind terminal 
velocity (v„) of ~350 km s"' , which is abnormally slow for its spectral type. We also analyse here more than 9 years of available 
RXTE/ ASM data. We study the long-term X-ray variability of the source and find it to be similar to that observed in the wind-fed 
supergiant system Vela X-I, reinforcing the idea that 4U 2206+54 is also a wind-fed system. We find a quasi-period decreasing 
from ~270 to ~130 d, noticed in previous works but never studied in detail. We discuss possible scenarios for its origin and 
conclude that long-term quasi-periodic variations in the mass-loss rate of the primary are probably driving such variability in 
the measured X-ray flux. We obtain an improved orbital period of Poit, = 9.5591 + 0.0007 d with maximum X-ray flux at 
MJD 51856.6+0.1. Our study of the orbital X-ray variability in the context of wind accretion suggests a moderate eccentricity 
around 0.15 for this binary system. Moreover, the low value of Voo solves the long-standing problem of the relatively high X-ray 
luminosity for the unevolved nature of the donor, BD +53°2790, which is probably an 09.5 V star. We note that changes in Vc„ 
and/or the mass-loss rate of the primary alone cannot explain the different patterns displayed by the orbital X-ray variability. 
We finally emphasize that 4U 2206+54, together with LS 5039, could be part of a new population of wind-fed HMXBs with 
main sequence donors, the natural progenitors of supergiant X-ray binaries. 

Key words. X-rays: binaries - X-rays: individual: 4U 2206+54 - stars: individual: BD +53°2790 - stars: winds, outflows - 
stars: variables: general - stars: emission-line. Be 



1. Introduction 

High Mass X-ray Binaries (HMXBs) are X-ray sources com- 
posed of an early-type massive star and an accreting com- 
pact object, either a neutron star (NS) or a black hole (BH). 
Depending on the n ature of the co mpanion, HMXBs are tradi- 
tionally divided ("see lCorbe3ll986l) into three groups: Classical 
Massive X-ray binaries, Supergiant X-ray binaries (SXBs) and 
Be/X-ray binaries (BeXBs). 

Classical Massive X-ray binaries are a very small group (in 
the Galaxy, only Cen X-3) of very bright (Lx ~ 10^** erg s^^) 
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persistent X-ray sources jLewin et alJll995l) . They have close 
orbits, with short orbital periods, and accretion is believed to 
occur through localised Roche-lobe overflow leading to the for- 
mation of an accretion disc. Orbits have circularised and, when 
the compact object is an NS, it has been spun up to short spin 
periods, as a result of angular momentum transfer from the ac- 
creted matter on to the compact object. 

In SXBs, the X-ray source is believed to be fed by 
direct accretion from the relatively dense wind of an 
OB supergiant, with little angular momentum transfer (see 
IWaters & van Kerkwiik,1989) . resulting in moderately high X- 
ray luminosities, Lx ~ 10^^ erg s ' . The orbital periods of these 
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sys tems are typically Pnrh ^ 15 d. Some orbits are almost circu- 
lar dCorbet & Mukaill20Q2h . but finite eccentricities have been 
measured for several systems, some of them rather high (e.g., 
e = 0.17 for 4U 1538-52; Clark 2000). The spin periods of 
systems with NSs are rather long, Pspin > 100 s. 

In BeXRBs, an NS orbits an unevolved OB star surrounded 
by a dense equatorial disc. These systems can be bright X- 
ray transients or persistent low luminosity Lx ~ 10^"* erg s"' 
sources. With the exception of a few peculiar cases like the 
microquasar LS I -h61 303 (Massi et al. 2004, and references 
therein) all of them appear to be X-ray pulsars. Their orbital 
eccentricities range from close to zero to very high, leading to 
the hypothesis that different kinds of supernova explosion are 
possible dPfahl et aljl2002al) . and their orbital periods are gen- 
erally of the order of some tens of days. Spin periods also cover 
a wide range an d there is a str ong statistical correlation between 
f orb and f spin llCorbellll986h . suggesting effective transfer of 
angular momentum from the material accreted. 

The peculiar HMXB 4U 2206+54 is difficult to place 
within this picture. It is a persistent source, with Lx - 10^^- 
10^^ erg s"' and variability on timescales of hours similar 
to those of wind accreting systems. However, it is one of 
the few HMXBs not displaying X-ray pulsations, although 
there is strong evidence that the compact object is an NS 
(Torreion et al. 2004; Blav et al. 2005). The X- ray flux is mod- 
ulated with a periodicity of 9.568 + 0.004 d dCorbet & Peeld 
l200ll) . which can only be interpreted as the orbital period. 
The mass donor in 4U 2206+54, BD +53°2790, is neither 
a supergiant nor a Be star, but a peculiar late O-type star, 
whose spectrum does not admit a standard spectral classifica- 
tion (Negueruela & Reig 2001). While most criteria favour an 
09.5 V star, there are some indications of a much heavier mass 
loss than expected for that spectral type, such as strong Ha 
emission, a P-Cygni profile in He ii /1 4686 and a very strong P- 
Cygni profile in the ultravio let C iv resonan ce doubl et, suggest- 
ing a more luminous star ( Negueruela & Reig-2001.;l'Blav et alJ 
I2OO6I) . If it has the luminosity of a normal 09.5 V star, it is 
located at a distance of ^2.6 kpc (if it were an 09.5 III star it 
would be located at ^4.8 kpc). 

In a recent paper we used INTEGRAL and VLA data to 
constrain the nat ure of the compact object in 4U 2206+54 
( iBlav et alJbOol . Based on existing radio/X-ray correlations 
for black holes in the low/hard state, we excluded the black 
hole scenario. On the other hand, INTEGRAL is the third X-ray 
satellite providing marginal evidence for the presence of a cy- 
clotron absorption line, leading to an NS with B - 3.6 x 10'^ G. 
However, two problems persisted in the NS scenario: the lack 
of X-ray pulsations and an X-ray luminosity one to two orders 
of magnitude higher than expected. 

Here we analyse an ultraviolet spectrum obtained with the 
International Ultraviolet Explorer (lUE), to better constrain the 
wind properties of the donor, and more than 9 years of data 
from the Rossi X-ray Timing Explorer (RXTE) to derive infor- 
mation on the mass loss from the optical star and the binary 
parameters. This work is organised as follows: in Sect. |2lwe 
analyse and model the lUE spectrum, in Sect.|3lwe present the 
RXTE data, in Sect.l^we study and discuss the long-term X-ray 
variability, in Sect.|5lwe focus on the orbital X-ray variability. 



and discuss the long-term wind variability in Sect.|6l We stress 
the existence of a population of wind-fed HMXBs with main 
sequence donors in Sect. and summarise our conclusions in 
Sect.il 



2. A measure of the wind terminal velocity 

The observed X-ra y luminosity of 4U 2206+54 is in the range 
-lO^^-lO^*" erg s~' ( iBlav et alJ2005l) . In contrast, the expected 
Bondi-Hoyle accretion luminosity for a canonical NS in a 9.6 d 
orbit around a low-luminosity 09.5 III-V star, with a typi- 
cally fast wind of ~1500 km s"', is of the order of or below 
10^"* erg s This value is critically influenced by the wind 
terminal velocity, Vco . To obtain a measure of Vk, for the wind 
of BD +53°2790 we have analysed the only publicly available 
high-resolution UV spectrum of this star, obtained with lUE 
on 1990 June 18-19, with a total exposure time of 20 ks (mid- 
dle time at JD 2448 061.59). This is flie high-dispersion lUE 
spectrum SWP 39 11 2, described in lNegueruela & Reid l EoOll) . 
but we have used the new reduction available at the INES^ 
database. A heliocentric velocity correction of 16.37 km s~' 
has already been applied to the source spectrum, and we have 
further applied a correction of - 62.7 km s~' to accou nt for the 
radial velocity of BD +53°2790 ('Abt & Bautz"'l963). We note 
that this value is the average of three measurements spanning 
55.0-72.2 km s"', obtained from September 1961 to May 1962 
and covering different orbital phases (the precise phase of each 
data point is uncertain due to the error in the orbital period 
and the huge timespan between these observations and the cur- 
rent ephemeris). Moreover, a preliminary radial velocity curve 
( Blav .2006.) shows total relative variations up to +30 km s~\ 
with a relative mean value of ~0+10 km s ' around the or- 



bital phase when the lUE spectrum was obtained. Therefore, 
we estimate that the radial velocity correction is accurate to 
-lOkms-i. 

We have followed the SEI (Sobolev with Exa ct Integration 
of the transfer equation) method, as outlined in iLamers et alJ 
( 1987) in order to calculate theoretical wind line profiles. This 
method is a modification of the So bolev approximation to the 
stellar wind problem (see ISobolevi.196 0). The procedure fol- 
lowed was to create a grid of theoretical profiles with differ- 
ent wind and photospheric parameters and then match them all 
against the observed profile. The indications from t he ext ensive 
study of the stellar winds of iGroenewegen et alJ ( Il989l) were 
followed and the parameters listed there for an 08 V star were 
chosen as initial values to build our grid of models. Theoretical 
models could not reproduce the data when using terminal ve- 
locities above 500 km s '. Terminal velocities in the range 
300-350 km s"' and a turbulent motion with a mean veloc- 
ity in the range 20-100 km s ' yielded line profiles which 
resulted in the best match against the observed profiles. We 
show in Fig.[l]an example of two UV lines from BD +53°2790 
matched against two theoretical profiles, calculated for 300 and 
350 km s-i. For the Civ AA 1548.19-1550.76 doublet an up- 
per limit for the turbulent velocity of 80 km s ' was found, 
while for the Niv /1 1718.55 line this limit was found to be 
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Fig. 1. Profiles of the observed UV doublet C iv AA 1548. 19-1550.76 {left panel, where v = km s ' coiTesponds to the blue line 
of the doublet) and the line Niv /i 1718.55 (right panel) with the theoretical models for 300 and 350 km s"' superimposed. The 
arrow in the left plot indicates Vbiack=350 km s"' . The wind terminal velocity is clearly slower than 500 km s ' . 



100 km s"'. We note that, due to the presence of many wind 
features in this spectral region and the relatively poor signal- 
to-noise ratio of the lUE spectrum, the position of the photo- 
spheric continuum is quite uncertain, and the choice of normal- 
isation parameters ca n affect our r esults. However, the C iv dou- 
blet is saturated, and lPrinia et al.i ( ■1990) showed that the point 
where saturated wind lines turn upwards towards the contin- 
uum level, namely Vbiack, is a good estimator of Vco- As shown in 
Fig.n]left, the C iv doublet provides Vbiack = Vco - 350 km s"', 
in good agreement with the above estimates by using the SEI 
method. Fina lly, we used the auto matic fitting procedure devel- 
oped bv Georgiev & Hernandez I tt2005 ), based on genetic algo- 
rithms, which provides Vco - 450 km s ' for the C iv doublet 
and Voo - 350 km s ' for the Niv line, also in good agreement 
with the values quoted above. 

Although an accurate measurement of the wind termi- 
nal velocity in BD +53°2790 is prevented by the relatively 
noisy WE spectrum, we stress that only models with Vtx, < 
500 km s ' could reproduce the wind line positions and 
widths. This value is much lower than the range of 1120- 
1925 km s^' for 09V stars o r the 1275-1990 km s^' range 
for 09.5 III stars jPrinia et"al][l99 0i). Therefore, we conclude 
that the stellar wind of BD +53°2790 is abnormally slow for 
its spectral type. A similarly slow wind of 400 km s"' has 
been measured in the 09.5 V extreme fast rotator HD 93521 
jPrinia etalJll99d excluded from their mean because it had 
very pecuha r profiles; see also iHowarth & Reidl 1 19931 and 
lMassalll995l) . Anothe r O-type star with an abnormally slow 
wind, of 510 km s-i jPrinia et alJll990l) . is HD 37022 (also 
kno wn as Ori C) , a spectral varia ble iri the range 04- 
7 V (IWalbornl[l98ll Emith & Fullertoiil EoOSl) . whose anoma- 
lous properties are interpr eted in terms of a misaligned mag- 
netic rotator. Interestingly, iBlav et alJ l EoOfil) have found that 
BD +53°2790 has a high rotational velocity, vsin; = 315 + 



70 km s , and have suggested a possible connection with 
HD 37022 based on their optical characteristics. 

3. The RXTEI ASM data 

In order to study the long-term X-ray behaviour of 4U 2206-1-54 
we have analysed X-ray data in the energy range 1.3-12.1 keV 
obtained with the All Sky Monitor (ASM) on board RXTE. The 
RXTE/ ASM data used here spans from 1996 late February to 
2005 late July (from MJD 50135 to MJD 53576), amounting 
to a total of 3441 day s or 9.42 years (the previous analysis by 
ICorbet & Pee'i3l2001 was performed with ~5.5 years of data). 
Each data point in the original Ughtcurve represents the fitted 
source flux of a 90 s pointing or 'dwell' on the source, with 
a mean of ^18.3 dwells per day in the case of 4U 2206-1-54. 
We have also analysed the one- day average lightcu rve of indi- 
vidual RXTE/ ASM dwells (see iLevine et alJI 1 996l for details). 
This lightcurve contains 3315 flux measurements, with data 
lacking only for 126 days (less than 4% of the total). Since 
4U 2206-H54 is a weak X-ray source, when there are only a 
few individual dwells per day the one-day average flux derived 
is not very reliable. To avoid spurious points with large error 
bars, we have also constructed a one-day average lightcurve 
with at least 5 dwells per day, which contains 3050 data points 
(92% of the one-day average data), and another lightcurve with 
at least 10 dwells per day, which contains 2608 data points 
(79% of the one-day average data). We will consider all these 
four lightcurves when performing the timing analysis, and we 
will refer to fliem as DBD (Data By Dwell), ODA (One-Day 
Averages), 5D-ODA (5 Dwell ODA) and lOD-ODA (10 Dwell 
ODA). 

The resulting lOD-ODA lightcurve is shown in Fig.|2 The 
source is clearly detected during most of the 9.42 year cov- 
erage, with a mean count rate of 0.36;^y 3*j count s"' (using 
weights as l/o"^; the plus and minus standard deviation of the 
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Fig. 2. RXTE/ ASM one-day average lightcurve of 4U 2206+54 after removing all the averages containing less than 10 individual 
dwells (lOD-ODA). Each panel represents approximately 3.3 years of data. Long-term X-ray flux variations on timescales of 
hundreds of days are seen. 



mean have been computed separately for the points above and 
below it). Assum ing a distance of 2.6 kpc to 4U 2206+54 
jBlav et alJ l2006l) . and taking into account that the average 
Crab count rate in the RXTE/ ASM is 75.5 count s we ob- 
tain a weighted mean and standard deviation of the absorbed 
X-ray luminosity of 3-12.1 keV) — 

(1.4+];^) (d/2.6 kpc)2 X 
10^^ erg s"'. We note that usin g a hydrogen column density 
of A^H = 11 X 10^^ ato ms cm~^ jTorreion et alJl2004 and the 
formalism described in Gallo et al. 72003"), the unabsorbed lu- 
minosity would only be around 5% higher. 



4. Long-term X-ray variability 

Long-term variability of the mean X-ray flux is clearly seen in 
the RXTE/ ASM data, as pointed out by Corbet & Peele ( 2001). 
To display this variability more clearly, we plot in Fig.|3^ the 
same data as in Fig. |2lbut after averaging (using weights as 
1/cr^) all data points within a running window of 30-day length 
(corresponding to ~3 times the orbital period). The count rate 
varies between 0.009 and 0.9 count s"', with a mean of 0.31 
and a standard deviation of 0.17. For comparison, the mean 
obtained with adjacent 30 d windows is 0.31, with a standard 
deviation of 0.16 and an error of the mean of 0.015, imply- 
ing that the variability is real at an 11-cr significance. Similar 
results are obtained by using running windows of 20 and 10 
days, although spurious points reaching even negative flux val- 
ues appear, due to the presence of intervals with few data points 
in the lightcurve. Therefore, we used the data shown in Fig.|3t 
to compute the range in luminosities, which turns out to be 
Ai.3-12.1 keV) ^ (0.035-3.5) (c//2.6kpc)2 x 10^^ erg s"'. This 



confirms the flux v ariations with a facto r of ~ 100 on timescales 
of years noted by iMasetti et alJ ( l2004l) when comparing data 
from EXOSAT, RXTE and BeppoSAX, but now obtained with 
the same satellite and detector (although the value of the mini- 
mum flux is uncertain in this case). 

This kind of lightcurve does not show any resemblance 
to that of any known Be/X-ray binary, persistent or tran- 
sient. It is, however, similar to those of wind-accreting X- 
ray binaries. As a comparison, we have plotted in Fig. ^jp the 
RXTE/ ASM lightcurve of the wind accretor Vela X-1, averaged 
and smoothed in exactly the same way as that of 4U 2206+54 
(after excluding data taken during X-ray eclipses, which other- 
wise translates into a mean reduction in flux of ~ 1.2 count s"'). 
Vela X-1 has rather similar orbital parameters to 4U 2206+54: 
it contains an NS (Pspin = 283 s) in a low eccentricity orbit 
(e = 0.0898, Pn rh ^ 8.9644 d) around the B0.5Ib star HP 77581 
(see Ouaintrell et aljE003l and references therein). As can be 
seen, the lightcurve obtained for Vela X-1 shows less variabil- 
ity (a factor ~4) but, like 4U 2206+54, it also experiences re- 
current peaks and troughs, with no clear periodicity (see next 
subsection). 

As noted bv lNegueruela & Reisl (1200 ll) and bv 'B lav et'aP 

(I2OO6'), BD +53°2790 is not a Be star, but a peculiar 09.5 V star 
with a relatively strong stellar wind. Therefore the comparison 
to the SXBs is fully justified, as the mass-loss mechanism in 
the mass donors is likely to be ide ntical, i.e., a radiative stellar 
wind (see Kudritzki & P ulsl200(i and references therein). The 
main difference between the two kinds of systems is the higher 
luminosity class of the mass donors in SXBs. This higher lu- 
minosity class will result in a higher X-ray luminosity through 
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Fig. 3. a) The data shown in Fig. |2l smoothed with a running 
window of 30 days of length, allowing a much better visuali- 
sation of the long-term flux variability, b) The same procedure 
applied to the RXTE/ ASM data of the wind-fed HMXB Vela X- 
1 (after excluding the data taken during its X-ray eclipses). 
RecuiTent peaks and troughs are seen in both lightcurves. 

two effects: 1) a higher mass-loss rate at the base of the wind 
and 2) a larger radius for the mass donor, which will place the 
NS closer to the surface of the OB star and hence in the region 
where the wind is denser and slower (although we have shown 
in Sect.lSthat the wind of BD +53°2790 is abnormally slow). 
Apart from this difference, all available evidence supports the 
idea that 4U 2206-1-54 is a wind accretor, not fundamentally 
different from the SXBs. 

We have inspected if there is any trend of hardening or soft- 
ening of the spectrum as the flux increases on these long-term 
timescales. We have done so by computing all possible hard- 
ness ratios that can be constructed with the three energy bands 
ofRXTE/ASM data (1.3-3.0, 3.0-5.0, and 5.0-12.1 keV). We 
have worked with the lOD-ODA lightcurve before and after 
smoothing. All the results obtained are compatible, within er- 
rors, with a constant hardness as the flux increases, although 
we emphasize that the poor statistics are a strong limitation of 
the analysis of RXTE/ ASM data of 4U 2206-H54. 

4.1. Timing analysis 

We have searched for periodic signals in all lightcurves using 
standard tec hniques l ike the Phase Dispersion Minimisation 
(PPM. iStellingwerd 1197 8) and the CLEAN algorithm 
llRoberts et al The periodograms obtained for the lOD- 

ODA lightcurve between 2 and 1000 d are shown in Fig. 0] 
As can be seen, the orbital period of ~9.56 d is clearly de- 
tected with both methods (with two subharmonics in the case 
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Fig. 4. Periodograms of the lOD-ODA lightcurve of 
4U 2206+54 obtained by using the PDM (top) and CLEAN 
{bottom) algorithms. Apart from the ~9.6 d orbital period, a 
trial period around 270 d is clearly detected by both methods. 



of PDM). On the other hand, significant signal is detected si- 
multaneously with both methods around trial periods of ~133, 
notably -267, -488 and -800 d. PDM also detects a 1-year 
signal, with a harmonic and a subharmonic, which is prob- 
ably the result of our window function, since it is not de- 
tected by CLEAN. Apart from the orbital period, these trial 
periods are also clearly detected when analysing the smoothed 
data shown in Fig. [3^. We note that a peak at a frequency of 
^ 4x10"'' d"' (coiTesponding to a period of 25 d) was already 
presen t in the power spectrum presented by ICorbet & Peele 



(2001) in their figure 2, and also noticed by IMasetti et al 
( 2004) . who questioned if it cou ld be a superorbital periodic- 
ity (see IClarkson et al.lE003albl for recent discussions on the 
topic). 

To investigate further the long-term variability, we have 
split the lOD-ODA lightcurve in 2 equal data sets spanning 
1720 d each, and then re-applied the PDM and CLEAN algo- 
rithms. The PDM results show a period in the range -255- 
265 d in the first part of the lightcurve, with a harmonic and 
two broad subharmonics. For the second part of the lightcurve 
there are minima around 180 and 415 d, and no significant sig- 
nal around 260 d. The output of the CLEAN algorithm reveals 
similar differences: 256 d, its harmonic and less significant 
peaks for the first data set, and different peaks at 90, 180 d, and 
many less significant peaks in the second data set. All this indi- 
cates that we are not dealing with a periodic signal of -260 d, 
but with a quasi-periodic one, present in the first part of the 
lightcurve but not in the second one, and that the flux varies on 
timescales of hundreds of days. This quasi-period of -260 d 
can be easily seen as alternative local maxima and minima in 
Fig.|3t during the first -5 year s of data (similar to the data anal- 
ysed bv lCorbet & Peelel200ll) . We note that very similar results 
to the ones discussed above are obtained when analysing the 
DBD, ODA and 5D-ODA Ughtcurves. 
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Fig. 5. Wavelet amplitude map of the lOD-ODA lightcurve of 
4U 2206+54. A quasi-period decreases from ~270 to ~130 d 
during the timespan of the data. 

A better understanding of how this quasi-period changes 
with time can be achieved with a two-dimensional time-period 
m ethod. For this purpose , we have used the method presented 
by'Szatmarv et al." ( 1994). We compute W{f, t), which we will 
call the wavelet amplitude. This value will be high if the sig- 
nal contains a cycle frequency / at the time r, and low other- 
wise. We have considered trial periods from 20 to 800 d, with 
a resolution of 1 d, and a total of 1720 times of analysis have 
been taken into account, which corresponds to one point every 
2 d. We show in Fig. |5] the wavelet amplitude map obtained, 
where white coloured areas mean significant trial periods and 
black areas non-significant ones, while the two grey triangu- 
lar areas have not been explored due to the presence of severe 
border effects (see RibQ_eLal . 2 001 for a detailed discussion of 
this issue). One can see the variation of the long-term ~200 d 
quasi-period over time, as well as the signal around 500 d and 
marginally at around 800 d. 

In Fig. |6l we plot, for any considered time r in Modified 
Julian Date, the trial period that displays the maximum wavelet 
amplitude. A decreasing quasi-period is present in the data, 
and excluding very short trial periods found at the limits of 
the data set and eventual jumps to the ~500 d signal, it varies 
from ~270 d at the beginning of the observations to ~130 d at 
the end. Very similar results are obtained when performing the 
same kind of analysis but using the smoothed data shown in 
Fig.H. 

In contrast, a similar analysis of the Vela X-1 data with 
PDM, CLEAN and the wavelet based method reveals no such 
kind of long-term quasi-periodic variability in this supergiant 
binary system. 

4.2. Superorbital quasi-period or wind variability? 

Superorbital periods with Pjup in the range 30-240 d and 
Psup/Porb values in the range 5-22000 have been found in 
a group of around 15 X-ray binaries (see .Wiiers & Pringle. 
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Fig. 6. Maxima of the wavelet amplitude map shown in Fig.|5] 
versus time for the lOD-ODA lightcurve of 4U 2206+54. The 
decrease of the quasi-period, with superimposed jumps, is 
clearly seen. 



119991: iQgilvie & DubuslEooH and references therein). These 
superorbital periods have often been explained as a preces- 
sion of the accretion disc due to warping induced as a con- 
sequence of illumination from the central source. In the case 
of 4U 2206+54, if the quasi-periodic signal were a superor- 
bital period, the Psup/Poih ratio would decrease from ~28 to 
~14 during the time interval covered by the RXTE/ ASM data. 
The variability of this quasi-period appears much higher than 
in other systems. Interaction of modes could be invoked to ex- 
plain such variations, and even the jumps to higher values of 
the period (Ogilvie & Dubus 2001). However, as discussed in 
■Torreion et al. (2004) , there is no evidence for the existence 
of an accretion disc in 4U 2206+54, and conversely, there are 
strong reasons to think that the X-ray emission originates be- 
cause of direct accretion on to the surface of an NS from the 
stellar wind, as discussed above and in Sect. |5l Therefore, we 
are inclined to think that the long-term X-ray variability we see 
is due, as suggested bv Masetti et aL (.2004) . to variations in the 
wind of the primary. 

EiTatic variability on timescales of hours has been detected 
with RXTE and BeppoSAX, and attributed to wind den sity 
inhomogeneities. While in the first case ^egueruela & Reid 
ll200lt) reported an increase of the 5-10/2-5 keV hardness ra- 
tio with incr e asing 2.5-30 keV intensity, in the second one 
Mase tti et alJ ll2004h report a constant behaviour for the same 
hardness ratio with increasing 2-10 keV intensity. This later 
result suggests that variations in the accretion rate due to vari- 
ations in the wind density do not change the slope of the soft 
X-ray spectrum. We have used the same scheme of hardness 
ratio versus intensity as in the later case with the RXTE/ ASM 
and found the same behaviour. Despite the large eiTors present 
in these data, this would be consistent with the hypothesis that 
long-term X-ray variability is due to a changing wind of the 
primary. 
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Interestingly, correlated variations between a variable 
emission component in Ha and X-ray flux have been ob- 
served on timescales of years in the wind-fed HMXB system 
LS 5039/RX J 1826.2- 145 (Reig et al. 2003; McSwain et al. 
I2OO4I: iBosch-Ramon et alJ (20051) . By analogy, the long-term 
variations in the X-ray flux from 4U 2206-1-54 might be due 
to changes in the wind of the mass donor BD +53°2790. 
Stochastic variability in the wind may easily explain variations 
in the X-ray flux on timescales of hundreds of days with no 
periodicity. It is, however, much less obvious what physical 
mechanism could result in the quasi-periodic variability ob- 
served in 4U 2206-H54. 

BD -t-53°2790 is too early to fit within any known cate- 
gory of pulsating stars - no 13 Cep star is known with a spec- 
tral ty pe earlier than BO d Lesh & Aiz enman^'l973*: 'Tian et al.l 
I2OO3I) . Moreover, the timescales of variability are much longer 
than in any B-type pulsator of any kind jWaelkens et alJl998l) . 
The possibility that non-radial pulsation activity is present in 
the Oef star BD -i-60°2522 has been sug ge sted to explain vari- 
ability on time-scales of hours (Rauwetal- 2003), but the vari- 
ability in BD -i-53°2790 has a timescale more typical of Mira 
variables. 

Although clearly BD -i-53°2790 cannot share a physical 
mechanism with Mira variables, there must be a physical rea- 
son driving the long-term quasi-periodic variability that we 
have detected. Long-term (quasi-)periodicities have been ob- 
served in a few other peculiar O-type stars, such as HD 108 
llNaze et alJ Ebol) . The recent discovery of a 538-d recur- 
rence in the spectra l changes of the Of?p star HD 191612 
llWalborn et alJ200 4') represents another interesting example of 
long-term periodicity of unclear origin. While for HD 191612, 
the possibility that a binary companion in a very eccentric or- 
bit drives the changes is tenable, the quasi-periodic changes in 
BD -t-53°2790 cannot be explained by a hypothetical third body 
in the system. 

Whatever the physical driver, if the changes in the X- 
ray flux are associated with variability in the mass loss of 
BD -h53°2790, the average X-ray flux of 4U 2206+54 may be a 
good tracer of the recent wind history of the donor. Taking into 
account the observed correlation between Ha emission and X- 
ray flux in LS 5039, it may be worth considering the possibil- 
ity that, at least in well detached HMXBs with wind accretion. 
X-ray monitoring might provide information on the long-term 
evolution of wind characteristics. 

However, the existence of chaotic short-term variability 
on timescales from minutes to days in the X-ray flux from 
4U 2206+54 (Neguer uela & Reidl200lt iTorreion et af]l2004t 
iMasetti et all 2004) prevents the comparison between available 
optical spectroscopy and X-ray flux measurements taken even a 
few hours apart. Clearly, a detailed multiwavelength campaign 
has to be undertaken to be able to fully interpret the data in 
each energy domain. 

5. Orbital X-ray variability 

ICorbet et alJ feOOO) performed a preliminary analysis of ~4.5 
years of RXTE/ ASM data of 4U 2206+54 (after correction 
from a wrong position by ~0.5° in the survey catalogue). 
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Fig. 7. DBD Ughtcurve of 4U 2206+54 folded using 
Porb=9.5591 d and fo=MJD 51856.6. Error bars represent the 
error of the weighted mean in each of the 20 bin per period, 
and not the standard deviation. Two orbital periods are shown 
for clarity. A quasi-sinusoidal pattern is evident. 



They found a quasi-sinusoidal modulation with a period of 
9.570+0. 004 d, and sugge s ted it c ould be due to orbital motion. 
Later on, ICorbet & Peeld ( I2OOII) performed a detailed analy- 
sis of ~5.5 years of RXTE/ASM data and refined this value to 
9.568+0.004 d. The folded lightcurve appeared again quasi- 
sinusoidal, and they derived an epoch of maximum X-ray flux 
at MJD 51006.1+0.2. 

5.1. Timing analysis 

As already explained in Sect. 14.11 we have performed a de- 
tailed timing analysis of the currently available 9.42 years of 
RXTE/ASM data. The orbital period is significantly detected 
in all data sets and using both the PDM and CLEAN algo- 
rithms. Since the PDM is an epoch folding method and there 
is clearly a great deal of long-term variability, we have pre- 
ferred to trust the results obtained with CLEAN. We note that, 
since the folded lightcurve is quasi- sinusoidal, CLEAN is wel l 
suited for this analysis (see, e.g., lOtazu etaP l2002l l2004h . 
The ODA and 5D-ODA lightcurves contain spurious data, and 
since the results obtained with CLEAN may be influenced by 
them (because weights are not used), we have not considered 
these cases. For the DBD and lOD-ODA cases, the results ob- 
tained are 9.5589+0.0005 and 9.5593+0.0005 d, respectively. 
Fitting a cosine function to the DBD lightcurve provides a pe- 
riod of 9.5591 + 0.0007 d and an epoch of maximum X-ray 
flux at MJD 51856.6+0.1. This will be the ephemeris con- 
sidered hereafter We note, however, that splitting the DBD 
lightcurve in two halves provides cosine fits with periods of 
9.5694 + 0.0014 d and 9.5520 + 0.0026 d, respectively, with 
a difference of 0.017 d, much higher than the error of the fit 
quoted above. This is probably due to the superposition of long- 
term variability. 
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Table 1. Computed orbital variability of L^cc for NSs (with Mx=1.4 Mq and Rx-lO km) accreting from the wind of different 
donors in eccentric orbits. The values for 09.5 V and 09.5 III donors (4U 2206+54) have been computed fo r different wind 
termi nal velocities. As a comparison we also show the results for Vela X-1 considering an inclination of 80° dOuaintrell et alJ 
l2003h . The observed RXTE/ ASM luminosities are around 1/3 of < L^cc >■ 



Spectral type 


(Mo) 


(Re) 


Mop, 

(Mo yr-i) 


(km s-') 




(d) 


e 


-^acc max 

(erg s-') 


^acc min 

(erg s-') 


Ratio in L^c 


(erg s-') 


09.5V 


16.0 


7.3 


S.OxlO^** 


350 


0.8 


9.5591 


0.15 


5.53x10'"' 


2.95x10'^ 


1.87 


4.06x10'^ 










500 


0.8 






2.19x10''' 


1.12x10'^ 


1.96 


1.57x10'^ 










1000 


0.8 






0.22x10" 


0.10x10'^ 


2.10 


0.15x10'^ 


09.5 III 


20.8 


13.3 


2.0x10^' 


350 


0.8 


9.5591 


0.15 


33.4 xlO" 


17.9 xlO'^ 


1.86 


24.7 xlO'^ 










500 


0.8 






15.2 xlO" 


7.54x10'^ 


2.02 


10.8 xlO'5 










1000 


0.8 






1.85x10'"' 


0.79x10'"' 


2.34 


1.21x10'^ 


Vela X-1 (80°, 1.96 Mq) 


24.2 


28.0 


1.0x10^*^ 


1100 


0.8 


8.9644 


0.0898 


5.47x10'* 


2.52x10'*' 


2.17 


3.75x10'" 



We show in Fig.0the DBD lightcurve folded by using these 
ephemerides. We have weighted the data according to 1/cr^ in 
each one of the 20 bins per period used, where cr is the er- 
ror of each individual dwell. As can be seen, the modulation 
is quasi-sinusoidal, but showing a slightly slower rise and a 
faster decay. This behaviour is unaffected by small changes in 
the Pofb used or by small shifts in Jq to average the data in a 
different way, indicating that it might be real. Nevertheless we 
must be cautious about this issue, since the long-term variabil- 
ity can affect the folded lightcurve in a noticeable way. On the 
other hand, the local minimum around maximum (phases 0.95- 
1 .00) is washed out when using slightly different values for Poib 
and/or to, which makes it more uncertain. 

As a further check of the orbital variability, and in an at- 
tempt to avoid the influence of the long-term variability, we 
have analysed the data in the following way. First of all, we 
have only considered the DBD data points falling on intervals 
of time when the average flux shown in Fig.|3i was in the range 
0.25-0.55 count s"', to avoid intervals of poor statistics and in- 
tervals of too much activity of the source. This corresponds 
approximately to 57% of the original DBD data. After that, we 
subtracted the average flux from the DBD data, to remove the 
remaining long-term variability. Finally, we analysed the re- 
sulting data set as previously done. CLEAN provides an orbital 
period of 9.5587+0.0005 d, compatible with the one reported 
above. The folded lightcurve obtained by using this filtered data 
set and the quoted period is very similar to the one shown in 
Fig. (but with a mean count rate of 0). The local minimum 
around phases 0.95-1.00 is also present, although again with 
a low significance. Regarding the orbital period, we prefer the 
former value, 9.5591 + 0.0007 d, because it was obtained with 
the original dataset, which contained about 2 times the number 
of flux measurements used in the latter case. 

5.2. Wind accretion in an eccentric orbit? 

The minimum and maximum RXTE/ ASM count rates in Fig.0 
are 0.26 and 0.46 count s"', respectively, implying a ratio be- 
tween them of 1.8. Translating the observed count rates into 
absorbed X-ray luminosities, we obtain an orbital variability 



covering the range L(i. 3-12.1 keV) - (1.0-1.8) (t//2. 6 kpc)^ x 
10^^ erg s"'. This degree of variability is expected in a wind- 
accreting system with a low or moderate eccentricity. In order 
to explore whether the observed variability can provide con- 
straints on system parameters, we have used a beta-law with 
spherical symmetry to model the wind of the donor, computed 
the position and velocity of the compact object in an eccentric 
orbit around it, and obtained the luminosity due to accretion by 
using a Bon di-Hoyle accretion model (Bo ndi & Hovldfl944t 
Bondi'1952'). A detailed explanation of the method is given 
in Reig et aL (2QQ3). A comparison between the equation we 
have used to estimate the X-ray luminosity due to accretion. 
Lace = GMxM ^ccJR-v, and that of the commonly used formal- 
ism bv Lamers et al.l(ll976i) . where Lace = ^M-^ccC^, reveals that 
the efficiency factor for the conversion of accreted matter to 
X-ray flux is ^ = GMx/Rxc^ = 1.48(Mx/M0)/(/?x/km). 

Since this is a very simple accretion model, we used 
Vela X-1 to check its validity. Folding the RXTE/ ASM data of 
the source with its orbital period, we find maximum and min- 
imum count rates of 6.4 and $3.2 count s"' (the presence of 
the X-ray eclipse close to the minimum X-ray flux prevents 
an accurate estimate in this latter case), thus providing a lumi- 
nosity ratio ^2. 0. Assuming a distance to the source of 1 .9 kpc 
dSadakane et alAl985J . these count rates translate into absorbed 
X-ray luminosities of 1.3 x 10^*= and 5 0.7 x lO""' erg s wifli 
an average of < L(i.3-i2.i keV) >- 1-0 x 10^^ erg s~'. On the 
other hand, by using logL/Lo = 5.53 dSadakane et alJll985l) 
and the relationship by Howarth & Prinial ( Il989l) . we obtain a 
mass-loss rate of M ^ 1 x 10"* Mq yr"'. For the wind terminal 
velocity we use the value Voo - 1100 km s"', obtained from 
lUE spectra jPrinia et alJll990l) . Hereafter we will use /3 = 0.8 
for the exponent of the wind law. The orbital p eriod and the ec- 
centricity are those of lOuaintrell et alJ (l2003l) . while we have 
considered their intermediate case of !=80°, which provides 
Mopt = 24.2 Mq, Rapt = 28.0 Mq, and Mx = 1.96 Mq (with 
the adopted values, ^ ^ 0.3). The expected accretion luminos- 
ity values are quoted in the last row of Tabled As can be seen, 
we obtain a ratio in Lace of 2.2, compatible with the ^2.0 ra- 
tio obtained with the RXTE/ ASM data. On the other hand we 
find < L(i.3-i2.i keV) >-< i'aee > / 4, since not all the ac- 
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Fig. 8. Computed variability of the accretion luminosity from 
the compact object in 4U 2206+54, assuming a neutron star 
with Mx = 1-4 Mq and Rx - 10 km, using a spherically 
symmetric wind and a Bondi-Hoyle accretion model. The vari- 
abihty is due to the orbital motion in an eccentric orbit with 
e = 0.15 around an 09.5V donor with Voc - 350 km s 
All parameters are quoted in the first row of Tabled] A quasi- 
sinusoidal pattern, similar to the observed one, shown in Fig.0 
is obtained. The observed soft X-ray luminosity is ~ 1/3 the ac- 
cretion luminosity. Phase corresponds to periastron, and the 
maximum of accretion luminosity is displaced due to the low 
wind velocity. 



cretion luminosity is released in the s oft X-ray band. Indeed, 
from BeppoSAX data of Vela X-1 bv lOrlandini et alJ flQQS). 
we found that between 1/3-1/4 of the total X-ray luminosity 
between 2 and 100 keV would be emitted in the RXTE/ ASM 
energy range. Therefore, the accretion luminosities obtained 
through our model yield a ratio of maximum to minimum lu- 
minosity compatible with the RXTE/ ASM data, and absolute 
values a factor of 3 to 4 above them, in agreement with obser- 
vations. 

In the case of 4U 2206+54 we have used Mx -1.4 Mq and 
Rx - 10 km for the NS (which provides ^ ^ 0.2). We have run 
models with the parameters corresponding to both an 09.5 V 
star and an 09.5 III star In the first ca se, we cons idered Mopt = 
16.0 Mq and R^pt = 7.3 Rq (Martins et alJl2005l). with M„„t ^ 
3 X 10 ~^ Mq yr"^^ (by using lo g L/Lp - 4.65 from Martin s etaP 
I2OO5I and the relationship bv lHowarth & Prinia .1989). In the 
second case, and using the same references, we have Mopt = 
20.8 Mq and /?opt = 13.3 Rq, with Mopt = 2 x IO^^Mq yr"'. 
We have tried different eccentricities and found that we can ap- 
proximately reproduce the quasi-sinusoidal shape of the folded 
lightcurve shown in Fig. by using e = 0.15, as can be 
seen in Fig. |8] for the 09.5V star considered above and us- 
ing Voo = 350 km s ^ We note that in this last figure phase 
corresponds to periastron, and the maximum accretion lu- 
minosity takes place at phase ~0.2 because of the low wind 
terminal velocity (i.e., if this model were correct, periastron 



would be at phase ~0.8 in Fig. 0. The results of our sim- 
ulations using an eccentricity of 0.15, for both the main se- 
quence and giant cases, and different values for the wind ter- 
minal velocity are quoted in Table ^ As can be seen in the 
main sequence case for Voo - 350 km s"', we are able to re- 
produce the 1.8 maximum to minimum X-ray luminosity ra- 
tio and obtain an accretion luminosity ^3 times the observed 
one (of < i(i.3-i2.i keV) 1-4 X 10^^ erg s"'), similar to 
the factor of ~2 b etween the 4-150 to 4-12 keV luminosities 
(Blav et al.l l2005l) . However, for higher wind velocities we un- 
derestimate the average X-ray luminosity, consistent with the 
low-speed wind observed by lUE (see Sect. |2}. In the giant 
case the system would be at 4.8 kpc, and the average X-ray lu- 
minosity < 3-12.1 keV) >= 4.8 X lO-'^ erg s"'. Therefore we 
should expect average accretion luminosities between 15 and 
20 X 1 0-'^ erg s" \ which can be achieved with Vco - 400 km s" ' , 
again compatible with the lUE data. We caution, nevertheless, 
that the wind terminal velocity has been derived with UV data 
taken several years before the X-ray flux measurements form 
RXTE/ASM. 

It is important to note that, for the unevolved nature of 
BD +53°2790 and the orbital parameters of the system, the low 
value of the wind velocity allows 4U 2206+54 to have an un- 
expectedly high X-ray luminosity even if containing an NS of 
1 .4 Mq as the accreting compact object. This solves the main 
problem in accepting the NS hypothesis, the other one being 
the lack of X-ray pulses that can be explained by simple geo- 
metrical effects (Blay et al. 2005). Finally, from the absence of 
X-ray eclipses, and for e - 0.15, we derive an upper limit on 
the orbit inclination of 82.5° in the 09.5 V case, and 77.5° in 
the 09.5 III case. 



6. On the long-term wind variability 

There are two simple ways to explain the long-term changes 
in the X-ray flux from 4U 2206+54: changes in the mass- 
loss rate of the donor and/or changes in the wind terminal 
velocity. Assuming a constant mass-loss rate of Mopt = 3 x 
10"^ Mq yr"', wind terminal velocity changes between 200 and 
1 100 km s ' should be invoked to provide accretion luminosi- 
ties in the range 10^'*-10''^ erg s"', necessary to explain the 
observed long-term X-ray variability discussed in Sect. |4] of 
i(i.3-i2.i keV) ^ (0.035-3.5) (£//2.6kpc)2 x lO''^ erg s"'. These 
velocity changes would lead to different orbital variability pat- 
terns: for the lower values of Voo we would have a higher X-ray 
flux with a maximum peaking at a phase 0.28 after periastron, 
while for the higher values of Voo the lower X-ray flux would 
have its maximum around phase 0.06 after periastron. To check 
if this is the case, we have split the DBD data in three intervals 
according to whether the 30-day averages shown in Fig.|3t are 
in the following ranges: > 0.5 count s"' (24% of data), > 0.3 
and < 0.5 count s"' (38% of data), < 0.3 count s"' (38% of 
data). We have subsequently folded these data with the orbital 
period, as done in Fig. We show the results in Fig. |9l For 
higher count rates the maximum X-ray flux takes place earlier 
than for lower count rates, contrary to what would be expected 
if only changes in Vco drove the changes in X-ray luminosity. 
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Fig. 9. The same as Fig.0but splitting the data in 3 different 30- 
d-average flux levels: high, middle, and low. The middle case is 
similar to the total shown in Fig.0 with a broad maximum cen- 
tered at phase 0. When the average flux is high, the maximum 
X-ray flux happens at phase ~0.85, while for the low flux case 
the maximum is around phase 0. 10. This behaviour is contrary 
to what would be expected if the average X-ray flux increases 
because of a decrease in the wind terminal velocity. 



If only changes in the mass-loss rate of the primary are in- 
voked, then for Voo - 350 km s ' we would need dramatic long- 
term variations in the range Mopt = 7 x 10^^-7 x 10"'" Mq yr"' . 
Moreover, a change in the mass-loss rate does not modify the 
phase of the X-ray luminosity maximum. Despite the relatively 
noisy RXTE/ASM data of 4U 2206+54, we can say that the 
scenario is not so simple. 

A general explanation of the various regimes of wind accre- 
tionon to a magnetised neutron star can be found in Stella et al. 
l|l986l). In our case, the accretion radius r^^c - ^GMx/v^^i 
jWaters & van Kerkwiiklfl%9i) varies along the eccentric or- 
bit, and takes values in the range 2-3x10" cm for the Vk, = 
350 km s"' case (~1.5xlO" cm for Voo - 500 km s"', and 
~0.5xlO" cm for Vk, = 1000 km s"'). This is one order of 
magnitude smaller than the separation of both objects at peri- 
astron passage. On the othe r hand, by using a surface magnetic 
field of Bo = 3.6x 10'^ G (.Blav et al.l2 005> and the dioole mag- 
netic field formula (B(r) = Bo{Rx/r)^', Waters & van Kerkwiik 
Il989 ) we can compute the magnetic radius along the orbit, 
given by r^ag - (B^R^/^npv^^^y , which takes values in the 
range 0.6-l.OxlO'" cm for all wind values considered. Indeed, 
considering the cases quoted in Tableland all possible orbital 
phases, the accretion radius is always between 5 and 40 times 
larger than the magnetospheric radius. Finally, since we do not 
know the spin period of the NS, we cannot compute the coro- 
tation radius. However, if we want to avoid centrifugal inhibi- 
tion of accretion, also known as the propeller mechanism (see 
IStella et aljfl98a) . we must have r^,- = (GMxPl^J^n^y^ > 
'"mag- Therefore, to have direct wind accretion in 4U 2206+54, 
the spin period has to be longer than 470 s. We have com- 



puted the magnetospheric radius as a function of the orbital 
phase, and it is at its minimum just before periastron and shifted 
in phase ~0.3 before the maximum of accretion luminosity. 
Therefore, we suggest that the behaviour shown in Fig.|5]could 
be the result of enhanced accretion when the magnetospheric 
radius is minimum, during epochs of high mass-loss rates of 
the primary. 

In addition, geometric variations in the X-ray irradiation 
of the wind of the primary can lead to changes in the ioniza- 
tion states of the wind sp ecies that we see, res ulting in varia- 
tions of the measured Voo ( iHatchett & McCra^ri9 77). This ef- 
fect has been observed in Vela X-1 ('K aneretalJh993 ) and in 
4U 1700-37 (Iping et al. 2004). Moreover, for higher X-ray ii-- 
radiation one has a higher ionization and a slower wind, thus 
a higher accretion rate and X-ray luminosity. The eccentric- 
ity of 4U 2206+54 could lead to additional significant changes 
in the amount of X-ray irradiation of the primary wind along 
the orbit, and thus in Voo, providing a different orbital variabil- 
ity pattern from the one obtained with a simple Bondi-Hoyle 
model. There are certain epochs when the RXTE/ASM maxi- 
mum to minimum orbital flux ratio appears to be a factor of ~5. 
This can be seen in data with relatively good signal-to-noise 
ratios after averaging with a 2-day running wir idow. An exam- 
ple of this variability is shown in Figure 4 of ICorbet & Peelel 
(2001). If this behaviour were only due to accretion variabil- 
ity in an eccentric orbit, a value of e ^ 0.4 would be needed. 
However, these factor of ~5 flux variations are not seen in the 
average folded lightcurve. Even when considering the split data 
in Fig. |9j variations with a factor in the range 2-3 are seen. 
Clearly, a radial velocity curve is needed to constrain the ec- 
centricity of 4U 2206+54 and allow a better modeling of the 
data. 

7. HMXBs with main sequence donors 

4U 2206+54 has been known since the early days of X-ray 
astrophysics ( Giacconi et al. 1972), although a complete un- 
derstanding of its accreting properties has only been possible 
after analysing an lUE spectrum. The abnormally slow wind 
of the donor in this binary system, ~350 km s"\ results in a 
relatively high X-ray luminosity, lO^'^-lO-'*' erg s"', for its rel- 
atively wide orbit, ~9.6 d, which has allowed its detection in 
all X-ray surveys. Iri addi tion, the system is relatively nearby, 
2.6 kpc ( Blav et alJEooi . and does not suffer an extremely 
high absorption, A^h = 1.1 x 10^^ atoms cm"^ diTorreion et alJ 
2004). Similar systems containing main sequence donors with 
normal fast winds would have one to two orders of magnitude 
lower X-ray luminosities, preventing their detection in existing 
X-ray surveys. 

The only exception i s the nearby (2.5 kpc) X-ray bi- 
nary system LS 5039 (see ICasares et alJl20Q5[ and references 
therein). The donor is an 06.5 V((f)) star dClark et al.|l200ll) 
displaying a fast wind of 2440 km s"' (McSwain et al. 2004), 
but the s ystem is present in the RO SAT All Sk y Bright Source 
Catalog dvoges et al.lll999t lMot(^h et al.lll997l) . because of its 
close orbit, ~3.9 d (Casares et al. 2005) , and relatively low ab - 
sorption, A^h - 0.7 x 10^^ atoms cm"^ ( Martocchia et al .120051) . 
We note that the X-ray luminosity of this system, ^ 10^"^ erg s^' 
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jBosch-Ramon et alJ2005l) . is on ly a small fraction, a factor of 
1/80, of the accretion luminosity dCasares et alJ2005l) . making 
its detection even more difficult. 

In this context, sensitive pointed observations of the 
Galactic centre with Chandra revealed the existence of a 
population of ~1000 low-luminosity hard X-ray sources 
(Wang et al., .2002: Muno e t al. 2003). Although a substan- 
tial fraction of them could be wi nd-fed accreting N Ss in 
HMXBs with main sequence donors (•Pfahl et al.'"2002b'), their 
nature is still iinder d ebate (see Bandvopadhvav et al. 2005 and 
Lavcock et alJ2005l) . 

On the other hand, available all-sky (or Galactic-plane) 
X-ray surveys are not sensitive and/or hard enough to detect 
sources of this kind a few kpc away from us, where interstellar 
absorption probably plays a crucial role. Therefore, the exis- 
tence of a population of wind-fed HMXBs with main sequence 
donors (and f spin > 100 s in the case of NSs to avoid the pro- 
peller effect), which would be the natural progenitors of SXBs, 
could be unveiled if new sensitive all-sky surveys with ener- 
gies above ~5 keV are performed (like the planned ROSITA- 
and EXIST^ surveys). Their study could have fundamental con- 
sequences for models of evolution and population synthesis of 
binary systems. 

8. Conclusions 

After a study of an lUE spectrum of BD -h53°2790 and -9.4 
years of RXTE/ ASM data of 4U 2206+54, we conclude that: 

1 . The ultraviolet spectrum reveals that, in mid June 1990, the 
09.5 V star BD +53°2790 had a wind terminal velocity of 
~350 km s"', abnormally slow for its spectral type, but sim- 
ilar to the ones measured in two peculiar fast rotators. 

2. An improved orbital period of 9.5591 + 0.0007 d was 
obtained, to be compared with the slightl y longer period 
of 9.5 68±0.004 d reported previously by Cor bet & Peeld 
1I2OOI ) using the first -5.5 years of RXTE/ ASM data. 

3. Long-term X-ray flux variability on timescales of hundreds 
of days is present in the data, compatible with a quasi- 
period that decreases from ~270 to ~130 d during the 
course of the RXTE/ASM monitoring between February 
1996 and July 2005. We conclude that this reflects changes 
in the wind of the donor, whose quasi-periodic nature re- 
mains puzzling. 

4. Using a Bondi-Hoyle accretion model, a spherically sym- 
metric wind with Voo - 350 km s"' and an eccentric orbit 
with e ^ 0.15, we are able to reproduce quite well the av- 
erage X-ray luminosity variations of the source with the 
orbital period, as well as the absolute X-ray luminosity of 
the system, which no longer poses problems in accepting 
a (non-pulsating) NS as the compact object in this binary 
system. 

5. The different patterns of the orbital X-ray variability for 
different average X-ray fluxes indicates that the long-term 
X-ray variability cannot be explained by variations in Voo 

^ ,http : //w ww ■ rssd ■ esa ■ int/index . php?pro j ect=Rosita, 
^ |http : // exist . gsf c. nasa . gov/. 



and/or M^pi alone. Due to the orbital eccentricity we expect 
changes in the magnetospheric radius and perhaps changes 
in Voo caused by variations in the ionization states in the 
wind of the primary due to changes in the X-ray irradiation. 
These parameters could vary as well on long-term scales 
for different mass-loss rates of the primary. We suggest that 
they might play a role in explaining the different patterns of 
orbital X-ray variability. 

6. If the magnetic field of the NS is Bq = 3.6 x 10'^ G 
( Blav et al. 2005), its spin period has to be longer than 470 s 
to allow direct wind accretion. 

7. Long-term coordinated observations in the optical and X- 
rays could confirm the proposed variability of the wind, by 
a correlation between X-ray flux and Ha excess. 

8. Observations are underway to obtain the radial velocity 
curve of 4U 2206+54 and constrain the eccentricity of the 
system, which is needed to properly model the observed 
orbital X-ray variability of the source. 

9. The nearby X-ray binaries 4U 2206+54 and LS 5039 are 
the only two known wind-fed HMXBs with main sequence 
donors. We suggest that more sensitive and harder X-ray 
surveys than the available ones could unveil a new popula- 
tion of objects of this kind, which are the natural progeni- 
tors of supergiant X-ray binaries. 
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